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ABSTRACT 


There are numerous examples of research dealing with towed and 
tethered bodies. Besides the classical problem of kites and towed gliders, 
recent work has been done on towed decelerators for reentry bodies and towed 
underwater devices. 

Previous analytical treatment of the body-cable system has assumed a 
rigid body problem, where the cable effect is accounted for by some force 
condition at the attachment point. The present approach treats the system as 
being essentially a cable problem, with the body dynamics giving end conditions. 

The mathematical form of the first order problem is a nonhomogeneous 
initial-boundary value problem in the partial differential wave equation. All 
equations, including end and auxiliary conditions, are linear with constant coef- 
ficients. Further, these equations uncouple to give a "lateral" problem and a 

"longitudinal" problem as in first order airplane dynamics. The solution of 

either problem takes the form of a transcendental characteristic equation for 
the stability roots. These roots are extracted by using an electronic computer 
and a roots locus plot. 

A series of tests on a tethered wind tunnel model provided a comparison 
of the theory with experiment. The equilibrium properties of the system were 
found by force measurements and photographs for various wind speeds, but the 
stability derivatives were found entirely from theory. This, along with the fact 
that the model was flying near stall, gave rise to a certain amount of difference 
between the theoretical and experimental results. Nevertheless, within the 
estimated error limits, the comparison is good and provides a convincing argu- 
ment for the theory. 
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INTRODUCTION 


The objective of this research was to devise a method by which one 
could predict the dynamic stability of a cable-body system subjected to fluid 
forces. Examples of such cable-body systems are numerous. Besides the old 
examples of kites and towed gliders, current examples are towed decelerators 
for reentering spacecraft (Fig. 15), towed underwater devices (Fig. 16), and 
towed and tethered "kite balloons" (Fig. 17). 

Analysis of cable-body system dynamic stability was initially directed 
toward kites and towed gliders, and the best examples of this work are given by 
Glauert (Ref. 6), Stfhne (Ref. 1 ), and Bryant, et. al. (Ref. 1). Current analy- 
sis, however, has been directed toward a wider range of cable-body system 
applications. For instance, besides towed glider research by Maryniak (Ref. 11), 
towed nonlifting surveillance devices were studied by Ftkin and Mackworth 
(Ref. 5), Mettam (Ref. 12), and Hopkin (Ref. 8). Also, a stability theory for 
towed underwater devices was developed by Laitinen (Ref. 9); and MacNeal 
(Ref. 10) derived a stability criterion for towed hypersonic decelerators. 

The approaches taken toward the problem by most investigators 
has been to consider the situation to be solely a body problem, with 
the cable accounted for by some force condition at the attachment point. 
This force condition takes the form of cable "stability derivatives", 
and is derived from the assumption that the cable is in an instantaneous 
equilibrium configuration with respect to certain of its end conditions • 
That is, the shape and tension of the cable are assumed to be either 
functions of the body's position, (Ref. l) , or functions of the body's 
position and velocity (Ref. ll). This assumption gives a good physical 
model for a large class of cable -body problems, in particular, those 
for which the cable is light and the system oscillations are slow. 

These conditions are met by most towed gliders 



and towed and tethered kite balloons. Further, the physical model has the merit 
of yielding a first order problem composed of a system of second order ordinary 
differential equations. These, in turn, give a polynomial for the stability roots. 

The present approach considers the dynamics of the cable, that is, the 
system is treated as a cable problem, with the body giving end and auxiliary 
conditions. Mathematically this means that the body f s ordinary differential 
equations of motion provide end and auxiliary conditions to the cable* s partial 

differential equations of motion. The result of the first order problem as 

shown in Chapter 3 is transcendental equations for the stability roots. The 

author feels that this physical model embodies the essential characterization of 

the cable-body system because the assumptions as shown later relate 

only to details about the cable construction, body construction, and fluid loads, 

and not about the mechanical nature of the cable-body system. An extension 
on the present approach is also included in reference 17* 
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1. THE CABLE EQUATIONS OF MOTION 


1. 1 The Complete Cable Equations 

The physical model of the cable considered is subject to the following 
assumptions: 

a. The cable is uniform along its length, perfectly flexible, inextensible, has 
a round cross-section, and is totally immersed in a homogeneous fluid. 

b. The Reynolds number of the cable r s crossflow is subcritical. 

The dependent variables of the cable are the coordinates of a point on it, 
£,y, and £>, and the tension, T, at that point. The independent variables 
are the distance along the cable, s, and time, t. The forces acting on an 
element of the cable, ds, are now found. Considering first the tension force, 
T, note that 

T= |T|n and f ' = |T + |jds |(iT+|^ds), 

thus, 

t t ,0 

dT = T' - T = ^Tn)ds. 

Further, 


9y~ * 3^r*- 

n ~ ds b l + 8s 2 + 8s b 3’ 


so, one has 


d? ’ ^ T i )dSe 'l + ^< T S» dS? 2 ■'^ <T S )dsb 3- <*• 11 

The gravity force, dFg, on the element, ds, is given by 

dF g = - pg(Cab^ +S5b 3 )ds. (1.2) 

Similarly, the buoyancy force, dF , on the element is given by 

J3 

dF_, = B(CSb +Sab )ds. (1. 3) 

13 1 o 

Treating now the fluid dynamic force, dF^, on the element, ds, note first that 


3 



where v is the cable element' s velocity relative to the reference frame (ft 
(see Fig. 1). Note that v^ is the relative velocity of the fluid to the cable 
element. Considering the plane of F^ and v^ t and introducing the cable 
coefficients, and C^, one obtains the components of the fluid force, dF^, 
as defined in Fig. 2: 

dF = C p |v I v Rds (1.4) 

a a 1 r 1 r ' 

and 

dF b = C b pv r X (v r X n)Rds. (1. 5) 

Now, 


so, 


7 s - Ifo + ft ? 2 

— ~ — 8v— 8^, -* 

V r = < U °* “ sK " < U& * + iT ) b 3 


r ' 8t 7 • at 2 ' ‘ 8t‘ 

Further, as described in page 3. 9 of Hoerner (Ref. 7), 

C = C +K(1 - C Z a) 3l/Z 
a a 

o 


and 


where 


C b = K(1 - C 2 <?)CS, 


Ca = 


(v r • n) 


( 1 . 6 ) 


(1.7) 


( 1 . 8 ) 


(1.9) 


Thus, (1.8), (1.9), (1.6), and (1.7) into (1. 4) and (1. 5) give the fluid-dynamic 
forces on the element. 

Equating all of the forces to the acceleration of the element, one has 


dT -.dP s + dF B + dF a + dF b = +S b 3 ) 4. 


( 1 . 10 ) 


so, (1. 1) — (1. 8) into (1. 10) give the complete cable equations of motion. Taking 


components, one has, for the b^ direction. 
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? y=i (T lf> + < ,R < c a *Ks 3 ;)t(ixs-f) 2 + ( f) 2 

Ot 0 


2 ■ 
+(US5 +||> ] (UCS - If) - KS 2 aa?pK X [{(US - 1|) 

- <us5 + f>ii> x <us - §*» - (pg - B)C5. 

Similarly, one has, for the e^ direction, 




2 2 
+ (IBS + |f> ] If - KS 2 5Q5pR x [{(UCS - |f) 


+( |> t(UsS+ |i )} a +{(lES .| ) | 


St 5s 


_ (U£ & 

^St ds ( M’xJxt J - 


dt 3s' 3t 


And finally, one has, for the b direction, 

u 


at o 


o o ' 

Sv ^ 5£> ~ 3C> 2s- w 

+ <ff) +(USa +^f) ] (USa +^-) - PRKS oCa 




- <f>S - < usS + f>0 <ua5; + f» - ft - ** 


(l.ii) 


(1.12) 


( 1 . 13 ) 
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1, 2 The First Order Cable Equations 

In order to simplify equations (1. 10) -*-(1. 13), a small perturbation 
analysis is performed such as to give the first order form of these equations. 
Consider, now, the £ axis to be aligned through the end points of the cable' s 
equilibrium configuration (Fig. 3). Further, consider a perturbation from 
equilibrium such that 

£ = £ 0 (s) +§(s,t), y = y’(s,t) and £ = £ q (s) + £’(s,t), (1.14) 

where £ q (s) and £ q (s) are the equilibrium values and £' (s, t),y' (s, t), and 
£,'(s,t) are the perturbation values from equilibrium. Also, 

T = T o +r(s,t), (1.15) 

where T q is the mean value of the equilibrium tension and r(s,t) is the 
perturbation value from equilibrium. Note now the important assumption made 
that the equilibrium tension is constant and equal to T q over the cable' s length. 
This approximation is closely realized for a large range of cable-body systems, 
for example, towed reentry decelerators and towed and tethered bodies where 

(pg - K)L s: T q and C^RL £ T q . (1. 16) 

Experimental evidence (Refs. 1,5, 12 and Chapter 4) has indicated that unstable 

modes of many such systems if they are present occur within the cable 

length where conditions (1. 16) are still satisfied. 

Now, assume small perturbations from the equilibrium position such 

that 


0[e], 0[e]U and T (s,t) = 0[e]T Q where e « 1. (1. 17) 

Note that it follows that 


ML 

3s 


= 0[e ] and 


f-=°[ £ V 


Proof: 

Considering the special case of motion in the y, £ plane, one has 

w 2 

dy = ds and d£'= edy, so, d£=e ds, 
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which gives 


M =e 2 Also 

5s * Al • at ^T* 

Allowing for motion in the I direction gives the same results. At this point, 
the assumption is introduced that the cable has a shallow curvature such that 


51 

a ~r ~ ** 0, and 1. (1.18) 

o as as ' ' 

This assumption is, in fact, consistent with conditions (1. 16). 

Now, taking (1. 17) and (1. 6) into (1. 9), and taking this, along with (1. 17) 

and (1. 18), into the complete cable equations, and further, dropping terms of 
2 

0[e ] and higher, one obtains the first order cable equations of motion. These 
are 


2 2 

= T o ^-^ +pR[(C & +KS 3 a)USaCa - 2KUS 3 oCa]-|^ 
at as o *■ 

+ pRU 2 K[S 2 ^C 2 a(3 - Sa) - S 8 a(S 2 a - 2C 2 £)]^-, 


(1.19) 


t 1 Ul 


p — = T^— ^ - pRU[C o +K(S 3 £ +C 2 Ss 2 a)]|f + 


at 


o. 2 

ds 


- pRKU^cfcSf^! 

as 


( 1 . 20 ) 


and 



= T 


a 2 ?,* 


as 


2 " pRUt(C a 


3~ 2~ 2~ 2~ at,' 

+ KS a)(l + S a) +KS aCahjj" 


- pRKU 2 [S 2 aCb(3Sa + 1) - 2S 4 aCa + 2C 3 aS 2 a]^. 


As an aside, the first order equilibrium equations are 


( 1 . 21 ) 
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( 1 . 22 ) 


T 1 + pRU 2 Ca(C +KS 3 a) - pRKU^afca +Cq?(E - pg) = 0 

3s o 


and 


T q j - pRU 2 S«(C a +KS 3 a) - pRKU Z S 3 aC?a +SS(K - pg) = 0. (1. 23) 

3s o 

Within the context of assumptions (1. 18), equation (1. 23) is the defining equation 
for the first order cable shape. 

1. 3 The Nondimensional First Order Cable Equations 
Define the following factors: 


at ~ ( " 


2U, 


)D 


5 = , 4U v n 2 a , y ~ r ~ - £' , y' . , s 

» 9 — ( o ) D » &nd ^,y, t, s . 


3t 


(1.24) 


One can now obtain nondimensional forms of equations (1. 19), (1. 20), and (1. 21). 
In particular, define the factors: 


o T b 

*£ o T _ , T 

C = O’ J " W' a 4pL ’ 


5 £Rb and . TS £lL 


4pU L 

Further, use these to define 

3 3 ^ 

k =J[(C +KS a)S oba - 2KS otCa] 
j. a 

o 


and 


kg s JK[S 3 aC 3 «(3 - Sa) - £? a(S Z a - 2C 2 a)J. 
Equation (1. 19) now becomes 


3s 


Similarly, define 


k. s J(C +KS 3 o +KC 2 oS 2 a) 
b a 

o 


(1.25) 


(1.26) 


(1.27) 


(1.28) 


(1.29) 


and 
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I 


k ? s JKCaS 2 «. 
Equation (1. 20) becomes 


, 2 ~ 


3s 


And finally, define 


and 


k 3 = J[(C +KS 3 a)(l +S 2 a) +KS 2 aC 2 a] 
o 


k^ = JK[S 2 QG?(3Sff +1) - 2S^3Ca(S^a - C 2 ®)]. 


Thus equation (1.21) becomes 


D 2 ? - C 2 ^-| + k D4 + k = 0. 
^~2 3 43s 

3s 


(1.30) 


(1.31) 


(1.32) 


(1.33) 


(1.34) 
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2. THE BODY EQUATIONS OF MOTION 


2. 1 The Force and Moment Equations 

The physical model of the body is considered to be subject to the fol- 
lowing assumptions: 

a. The body is rigid, symmetric with respect to the n , n plane (see 

X o 

Fig. 5), and completely immersed in a homogeneous fluid. 

b. The cable is perfectly free to pivot at the attachment point. 

c. The center of buoyancy is on the n^ axis. 

It is important to note that the body equations will be expressed in terms 
of x,y,z, and the Eulerian angles <r, 0 and <p relative to (ft. Although this is 
different from airplane practice, these coordinates are necessary in order to 
relate the body equations to the cable equations. However, the force and moment 
equations are derived relative to n , n 0 and n because fluid-dynamic effects 
on a body are traditionally taken in these directions. But, by transformation 
equations, these force and moment terms are eventually expressed in terms of 
x,y,z,t/r , 0 and <p. 

The force-acceleration equations are, as in Etkin (Ref. 4), 



Fj = m(u +qw - rv), 

(2.1) 


Fg = m(v + ru - pw), 

(2.2) 

and 

F g = m(w +pv - qu). 

(2.3) 

And similarly, the moment-angular acceleration equations are 



M 1 = ’xx p - 'xz p t « I ra r - ! xz p > - rI yy q ' 

(2.4) 


M 2 = W 1 +r(I xx p - 'xz r > - « I zz r - ’xzrt- 

(2.5) 

and 

m 3 “ y - ! xz p +i y y pq ■ q(I xx p - y>- 

(2.6) 

Note that u,v, and w are defined by the velocity of the mass center: 
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V c = Un i +Vn 2 +W V (2.7) 

Also, p,q, and r are defined by the body' s angular velocity with respect to (R: 

«- pn 1+ qn 2+ riT 3 . (2.8) 

Consider now the relations between e. and n! based on the Eulerian 

1 1 

angles as defined in Fig. 6. These are 

e = CfCdn +(C4iS6S<p - StyC^rT + (C$S0C<t> +Sl/iS(f>)n , 

11 A o 


and 


e = SfCdn +(Sij!SeS(j) + CfC<p)n + (SijiS0Sip - C$S0)n , 

A 1 A o 


e = - S0n +C9S4>n + C9C(pr\ . 
o 1 A o 


These relations give, by virtue 


v c = xe l +ye 2 +ze 3 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 
( 2 . 12 ) 


and equation (2. 7), that 

u = xCijiCe +ySilfC0 - zS0, (2. 13) 

v = x(C'i|fS0S0 — SfCcp) +y(SjlrS0S(p + CfC<p ) +zC0S<p, (2.14) 

and 

w = k(CTllS9C<p + StfScp) +y(S0SilfC 4> - CiliS<p) + zC 9C(f>. (2. 15) 

Similarly, resolving (2. 8) with (2. 9), (2. 10), and (2. 11), one has 

p = "<p- $S0, (2. 16) 

q = 0C<p + $C0S<£, (2.17) 

and 

r = fC0C<p - 0S(j). (2.18) 

Equations (2. 13)-*(2. 18) into (2. l)-»(2. 6) give the force and moment equa- 
tions in terms of x,y, z,f, 0, <f>, and their derivatives. These equations are non- 
linear, but in the spirit of the stability analysis, a small perturbation is per- 
formed, and linear, first order equations are derived. 


2.2 The First Order Force and Moment Equations 

Consider a perturbation of the Eulerian angles and their derivatives such 

that 

11 
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(2.19) 


0 = 6 q + 0 , <t> = <P Q + <t>, if' = f Q + ft 

'e=d Q +0, ^>=^ o +<^, and f=f o +f, (2.20) 

where the "o" quantities are reference values, and the quantities are the 
perturbation values. Further, define the reference configuration of the body to 
be that of static equilibrium, thus, 


eh = ft = ft = e - b =o, 

o o o o o 

and 0 is a fixed value according to the condition that the n, axis passes 
o 1 

through the attachment point and the mass center (see Fig. 5). 

Now, assume small perturbations such that 


( 2 . 21 ) 


f ♦ • 

?,$,$ = 0[e], x, y, z = 0[e ]U, and = 0[e ](U/b), where £ « 1. (2.22) 

Substituting (2. 13)-*(2. 22) into (2. l)-*(2. 6), and dropping those terms containing 
2 

an 0[e ] or higher, one obtains the first order form of the force and moment 
equations. These are 


F, = m(xC0 - z'S0 ), 
1 ' o o 

F 2 = my, 


F n = m(zC0 +xS0 ), 
3 ' o o 


M 


= I % - (I S0 +1 C 0)|, 

1 XX XX o xz o 


M = I 0, 

2 yy 


and 


(2.23) 

(2.24) 

(2.25) 

(2.26) 
(2.27) 


M = (I C0 +1 80 ft -I 3. (2.28 

3 v zz o xz o y xz 

Finally, the dynamics of the body deals with its motion with respect to the inertial 
reference frame, dt. However, the fluid dynamics of the body depends on its 
motion relative to the fluid stream, (R 1 , which is 


(v ) = v - Ue, = u n +v n 0 +w n . (2.29) 

v c'r c 1 rl r2 r 3 

For motion subject to the small perturbation conditions, (2. 19)-»(2.22), equations 
(2. 13)-*(2. 15), and (2. 29) give 
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(2.30) 


u = xC0 - zSf) - UC0 , 
r o o o 

v r = tty +y - US0 q ?, (2.31) 

and 

w = zC0 4-xS0 - US9 - UC0 0. (2.32) 

r o o o o 

Now, considering again the equilibrium reference condition, one defines velocity 

perturbations by 

u = - UC0 + u’ , v = v' , and w = - US0 +w’ . (2. 33) 

r o r r o 

Thus, from (2. 30)-»(2. 32), (2.33) gives 

u' = kCQ - zS 0 , (2. 34) 

o o 

v’ = U# +y - US0 o <?, (2.35) 


and 


w' = zC0 4 xS0 - UC0 0. 
o o o 


(2.36) 


Similarly, the body's acceleration with respect to the fluid stream is given by 

- c - 

a = ~ — -fa + w x v , 
c dt c 

r 


which gives 


(ft 


a 

c 

r 


dv 

df a Vi + V2 + Vs- 


(2.37) 


And, for the small perturbation case, (2. 19)-*-(2. 22), equations (2. 13)-*-(2. 15) 
and (2. 37) give 


a, = xC0 - zS0 , 
loo 


a 0 = y. 


(2.38) 

(2.39) 


and 

a =zC0 4xS0 . (2.40) 

3 o o 

Note that the acceleration of the body with respect to the fluid stream, <R' , is 
identical to that of the body with respect to the inertial reference frame, (R. 
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Finally, one obtains the small perturbation forms of the angular accelera- 
tion components by taking (2. 19)-*(2. 22) into (2. 16)-*(2. 18). This gives 

p=^-^S6» o , (2.41) q=0, (2.42) and r = Cdjj). (2.43) 

2.3 The Nondimensional Form of the First Order Force and Moment Equations 
The factors used to nondimensionalize the terms in the force and moment 
equations are identical to those used in American airplane analysis, except that 
no distinction is made between a "longitudinal" and "lateral" characteristic length. 
Now define 


A a F .F >F >B,mg,T 

C x ’ C Y» C z ’ B ’ mg ’ T o H 2 . 

x Y z ° (pUS/2) 


(2.44) 


C„C ,C b 


M 1' M 2* M 3 


1 m n (pU 2 Sb/2) 


(2.45) 


^ A A A 


X ’ y ’ Z > R B’ R a A.. u-.V.W 
x,y,z,2R,2R a = ^ , u,v,w = — , 


(2.46) 


a i ,a 2 ,a 3 “ ,„„2 


a i» a 2 » a 3 * a a p,q,r 
, P , q»r = - ili 


(2U /b) 


(2U/b)’ 


(2.47) 


, 2U, 4m 

t _ -g-t> n = pgb , 


(2.48) 


I ,1 ,1 ,1 

xx’ yy> zz xz 

xx* yy’ zz ’ xz " ps (b /2) 3 ’ 


(2.49) 


D < (2-50) and D 2 ( 


Introducing these into the force and moment equations, (2. 23)-*(2. 28), gives 

2 a 2 a 

C = x - /iSO^D z, 

C v = nD y, 


(2.51) 

(2. 52) 
(2. 53) 
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(2.54) 


C Z = M S0 o d2 ^ +/iC6> Q D 2 z, 

c = i d 2 ? - (i so +i ce )d\ 

r. xX ' w n v7 n' ■ 1 


XX O XZ O 


and 


C =i D^, 
m yy 


C = (i C0 + i SO )D 2 2) - i D 2 $. 
n zz o xz o' Y xz 

Further, equations (2. 34)-*(2. 43) become 


(2.55) 
(2. 56) 


(2. 57) 


u = C0 Dx - S0 Dz, 
o o * 

(2.58) 

v = Dy + $ - S e Q $, 

(2.59) 

w = S0 Dx + C0 Dz - C0 0, 

O O 0 

(2.60) 

A 2 A 2 A 

a = C0 D x - S0 D z, 
1 o o 

(2.61) 

A 2 A 

a 2 = D y, 

(2.62) 

* 2^ 2* 
a = C0 D z + S0 D x, 
3 o o 

(2.63) 

p = D^> - S0 q D$, 

(2.64) 

A 

q = D0, 

(2.65) 

r = C0 q D$ 

(2.66) 


2.4 The Force and Moment Coefficient Terms 

Consider the forces, F., and the moments, M., (i = 1,2,3) on the 
body. Rewrite these as 

F. = F. + AF. and M. = M. + AM., (i = 1, 2, 3) 
o o 

where F. and M. are the reference values, and AF. and AM. are the 

ii ’ii 

perturbed quantities. Defining the reference condition to be that of static 
equilibrium, one obtains 

15 
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(2.67) 


F. = M. = 0. So, F. = AF. and M. = AM.. 

11 li li 

o o 

Now, in the spirit of small perturbations, AF. and AM. are assumed to vary 
linearly in the perturbed velocity, the acceleration, and the angular velocity of 
the body relative to . Also, accounting for the body weight and buoyancy, one 
assumes that AF. and AM. vary linearly with the perturbed Eulerian angles. 
Finally, a cable force and moment contribution is accounted for by AF and 
AM p ,. Thus, the general expressions for the force and moment terms are: 

9F. 9F. 9F. 9F. 9F. 

AF = — -u ' + — ~v ' + — ’ h -a h 

i 9u' 9v' 9w' 9a 1 9a 2 

J- 


9F. 9F. 9F. 9F. 9F. 

+ 9a a 3 + 9p 13 + 9q^ + 9r r + 9 

tj 


( 2 . 68 ) 


and 


9F. 9F. 

+ ^ + l^ +AF Ci < i=1 ’ 2 ’ 3 > 


9M. 9M. 9M. 9M. 9M. 

l t l f l l l 

~ 9u' U + 9v ,V + 9w ,W + 9a, a i + 9a„ a 2 


9M. 9M. 9M. 9M. 9M. 

+ 8^3 ^ + aT 1 + '5r r 


(2.69) 


9M. 9M. 

+ 8T ? + 1T' S * AM CI < i = 1 ’ 2>3) 

Further, when one defines the nondime ns ional parameters. 


„ 9F. 9F. 9F. 

C c c = <—)— — — 
Fi ’ Fi ’ Fi _ 'pUS ; 9u' ’ 9v’ ’ 9w’’ 
u v w 


(2. 70) 


„ 9M. 9M. 9M. 

C c c - (- 2 -) - - - 

Mi ’ Mi ’ Mi = 'pUbS' 9u’ ’ 9v’’9w'’ 

U V w 


(2.71) 
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„ 9F. 9F. 9F 

c C C = /— v - L L 

Fi. ’ Fi.’ Fi. “ VSb ) 9a, , aa 0 , 9a n ’ 
u v w 12 3 


, 9M. 9M 9M. 

C C C = ( — — — V - 

Mi. ’ Mi. • Mi. - ' ou 2'9a., ’ 9a„ ’ 9a 0 ’ 
u v w pSb 12 3 


. 9F. 9F. 9F. 

C c C = ( — — — \ - — i 

Fi * Fi ’ Fi ~ 'pUSb* 9p ’ 0q ’ 9r ’ 
P q r 


„ 9M. 9M. 9M. 

Q (-; o - t — ) i 1 

Stt '°Mi " t .™2 9p * 9q * 9r * 

p q r pUbb 


„ 9F. 9F. 9F. 

2.1 i l 


°Fi * °Fi ’ C Fi " ( 2* 9$ * 90 ’~9^’ 

>/) e <p pu s ^ 


„ 9M. 9M. 9M. 

C C C - /- 2 * 1 1 1 

Mi ’ Mi ’ Mi = ( 

0 0 0 pUSb 


2 )_ 9^’ 90 ’ 9? 


and 


2AF 2AF 

C Fi - 2 = CLR 

C pl^S pU 2 !^ S 


2AM 


Mi, 


2AM. 

l c l clr 

\ c 


C pU^Sb pl^LRb s 
equations (2. 68) and (2. 69) become 


C = r u C v + C w+C.a+C a+C a 
Fi Fi Fi Fi Fi 1 Fi 2 Fi 3 

U V w U V w 


* C Fi P +C Fi 5 +C Fi f +C Fi * +C Fi/ +C Fi/ +C Fi 
p q r 0 0 0 C 

(i = 1,2,3) 


and 


(2.72) 


(2.73) 


(2.74) 


(2.75) 


(2.76) 


(2. 77) 


(2. 78) 


(2.79) 


(2. 80) 
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(2. 81) 


C = C . u 4 C.. ... v + C_ . w+C a, +C„. a 4C . a„ 

Mi Mi Mi Mi Mi. 1 Mi. 2 Mi. 3 

u v w h V w 

+ C p +C. t , q + C . r +C $ +C 0 +C... $ + C... 

Mi Mi Mi Ml, Mi. Ml, Mi 
p q r ft 9 4> C 

(i= 1,2,3) 

Now, as in airplane stability analysis, cross derivative terms are dropped; 

that is, fluid dynamic stability derivatives of symmetrical quantities with respect 

to unsymmetrical variables are dropped. Thus, 

C ,C ,C ,C , C ,C ,C V , and C v = 0. (2.82) 

m m. m m X X. X X 

vvprvvp r 

Similarly, fluid dynamic derivatives of unsymmetrical quantities with respect to 

symmetrical variables are dropped. So, 


(2. 83) 


F +B = (B - mg)e 
g J 


C Y * C y > C y.’ C Y.’ C Y ,C Z ,C Z. ,C Z ,C Z ,C i ,C i ’ 
uv uwqvvpruw 

C „ , C „ , C „ , C , C , C , C , and C =0. 
i. 1. £ n n n. n. n 

uwquwuw q 

Further, considering the gravity and buoyancy effects, one has 

(see Fig. 7), 

which, by (2. 11), gives 

F + B= (B - mg)(-S0ri^ + COS <pn 2 + COC^n^). 

Further, the perturbation value from equilibrium is found by using (2. 19)-^(2. 22), 

2 

and dropping terms of order e and higher. This gives 

A(F 4 B) = (B - mg)(-C0 0n. 4 CO 3>iT 0 - SO ?nj. (2.84) 

g o 1 o 2 o 3 

Thus, (2.84) and (2.76) give 

V C X 'V c Y' c z; “ d c z -° (2 ' 85) 

Ip (p ip 0 Ip <p 


and 


C v = -(B - mg) CO , C v = (B - mg)C0 , and C = -(B - mg)S0 (2. 86) 

A ^ O i o Zj u 

6 $ 0 
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Similarly, the moment about the mass center due to buoyancy effects is given by 

M B = R B x B = r b B ^i X V (see Fi S* 7 >* 

Using (2. 11), one obtains 


M b = R B B(-C0C<^n 2 + C0S0 iT 3 ), 

for which, again, the perturbation value from equilibrium is found by using 

2 

(2. 19)-*(2. 22), and dropping terms of order e and higher. This gives 


= R B B(£ V n 2 + Cfl o?V* 


Thus, (2. 77) and (2. 87) give 

C , C , C , C , C . , C , , and C . = 0, 
m, m,’ n ’ n, i, £’ i, 

0 e 0 e 0 


(2. 87) 


(2. 88) 


C =RBSe and C = KBC0 . (2.89) 

m e on <p 

Now, the cable effects on the body are considered. The cable terms in the 
force and moment equations provide the mathematical link between the body' s 
motion and the cable' s motion. The cable force is 


-► 3x 8v -> 3z -* 

T a = - T o^ )e l +( ^ e 2 + W e 3l* 
a a a 


(2.90) 


Using equations (2. 9)-*(2. 11), one may resolve this into the n^n^, and n^ 
coordinate directions. Further, consider a perturbation of the cable from 
equilibrium such that 


.3x v .3x v 3x' v .3y. ,3y\ , / v4j \ / v<j v 

W =( al ) +( ai' ) ’ W = ( aP and te* = +( “ ) * 


(— ) 

W 


,9x' 

[ 9s 


'9s' 


(— ) 

w 


(-) 

l 9s ; 


' dS 


(2.91) 

aaaaa aaa 

where the ", " quantities are the equilibrium values, and the primed quantities 

are the perturbation variables. Consistent with the previous small perturbation 

analysis, (1. 17), the primed terms are said to be of order e. Thus, substituting 

2 

(2.91) into equation (2.90), and dropping terms of order e and higher, one 
obtains the first order cable force. Now. the equilibrium cable force is 
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(2. 92) 



a a a a 


and, subtracting this from the first order cable force gives one the expression 
for the perturbation value of the cable force from equilibrium. Nondimensionalizing 
by (2.78), one obtains the components of this expression: 


a 32' 9x' 

: v = t [(— ) so - (—) ce + 

o 3s o 3s o 

C a a 


{<S> S9 o + <§> Ce a> ?1 - 

a a 


(2.93) 


v ■ -v<tK° + (<f> se o 4 <S > ce ^ - <§> *>• 

C a a a a 


(2. 94) 


and 


^ 3z' 3x' r 3x 3z, ^ 

C Z * " T o l( “ei‘ ) Ce o 4 <1T> S9 o 4 fe> Ce o - « se >- < 2 ' 95) 

C a a a a 

In a similar fashion, the cable moment terms are derived. The moment 

on the body due to the cable force is 


M = T X R n - (2.96) 

C a o 1 

where T is given by (2. 90). Again, when one uses equations (2. 9)-*(2. 11), 
a 

equation (2. 96) may be resolved into and n^ components. And further, 

as by (2.91), a small perturbation from equilibrium is taken, and terms of order 
2 

e and higher are dropped. This gives the first order form of the cable moment. 
Now, the equilibrium cable moment is 


M = -R T [(|^) S3 + (|^) CO ]n 
C a o v 3s o v 3s' o J 


2 ’ 


(2.97) 


o a a 

and subtracting this from the first order cable moment gives one an expression 
for the perturbation value of the cable moment from equilibrium. Nondimen- 
sionalizing by (2.79), one finds that the components of this expression are 


C 


X 


c 


= 0 , 


(2.98) 
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(2.99) 


c m - - “o^faT 1 C9 o 4 {<f» ce „ - <H> S9 >- 

C o 


0Z* 


-9x , 


,9z v 


a 


a 


and 


C „i, * ft .*o«*K 4 {<S> S »o 4 <S>_ Ce o>5 - <§> a. 

V 


( 2 . 100 ) 

a ~ a a "a 

Finally, substituting (2. 82), (2.83), (2.85), (2.86), (2.88), (2.89), and 
(2. 58)-*(2. 66) into (2. 80) and (2. 81) gives one the force and moment expressions 
in terms of the x,y,z,$, 0, $ coordinates, namely, 

c v =■• [(C ce + C S 0 )D 2 + (C C0 +C S0 )D]x +[(C C 0 
X X. o X. o X o X o X. o 

u w u w w 


- C X. S6 o> d2 4 < C X Ce o - C X S V D > 2 + ‘ C X D - C9 o( C 3 


w 


u 


w 


+ (B - mg)}]3 + C , 


( 2 . 101 ) 


c y = (c y d +c y D)? +[(C 0 o c y - se o c Y )D + C y $ 

v v r p v 

+ [C y D - {S0 o C y - (B - mg )C(T}]$ + C y , 
p v C 


( 2 . 102 ) 


C r 7 = [(C ry C 0 rt - c 7 SflJD +(C Z C0 o ^C z s 6 o )D]Z 


Z. o Z. o' 
w u 


w 


u 


*l< c z , Ce o + C z. s V D 4 < c z ce 0 +c z S V D|X 

u w u w 

+ [C z D-{c z C 0 o +S 0 o (fe -mg)}]?4C z , 
q w C 


(2. 103) 


C i =<C *. D +c « D)y + (C D-C SO o )0 

V V p v 


+ [(C ce - c S0 )D + c 30, 
*r ° p v 


(2. 104) 
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C = UC CO - C SO )D +(C CO - C S6 )D]z 
m m . o m. o m o mo 

w u w u 


and 


+ [(C CO +C SO )D + (C C0 +C S<9 )D]£ 
m. o m. o m o m o 

u w u w 

+ [C d - (c ce 4bs e \]0+c 

m k m O O J m _ 

q w C 


C = (C D + C D)y + [(C C6 - C S6 )D + C $ 
n v n. n n o no n 

v v r p v 

+ [C D+<TrBC0 -C Sd}]$ +C . 

1 n V. o n oJ n _ 

p v C 


(2. 105) 


(2. 106) 


Thus, equations (2. 93)-*{2. 95) and (2. 98)-*(2. 100) into equations (2. 101)-*(2. 106) 
and these, in turn, into equations (2. 52)-*(2. 57) constitute the complete force 
and moment equations for the body. Note that the fluid dynamic force coefficient 
terms, (2. 70)-*<2. 75), may be directly related to the "stability derivatives" of 
standard airplane notation. The transformation equations to relate one to the 
other are given in the Appendix. 


2. 5 The End and Auxiliary Conditions as Given by the Force and Moment Equations 
As mentioned in the introduction, the key to the solution of the cable-body 
problem is to solve the cable equations, where the body equations of motion provide 

end and auxiliary conditions. To this purpose, the body equations of motion must 

be rearranged and combined so as to isolate the cable terms. First, note that 
(2. 94) and (2. 100) combine to give an auxiliary condition: 

C 4R n = 0 (2. 107) 

n c a Y c 

A second auxiliary condition is given by (2. 95) and (2.99): 

C - R C =0. (2. 10S) 

m c a z c 

Also, a third auxiliary condition is given by (2.98): 

C =0. (2. 98) 
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Further, (2.94) gives an end condition: 


'as' o 

a 


(2. 109) 


Similarly, (2. 93) x S8 q - (2. 95) x CO gives a second end condition: 


<1r> - r <c x_ S8 o - c z_ c9 o> - q 

a T 


.«) ? 
a 


( 2 . 110 ) 


And finally, a third end condition is given by (2.93) X C 0 q +(2.95) X S0 q : 


<T5T> + , % C9 o +C Z r S9 o» + <S> ? - 

a T C C a 

o 


( 2 . 111 ) 


Now, these conditions may be expanded into full form by using equations (2. 101) — 

(2. 106), (2. 93)— (2. 95), (2. 98)— (2. 100) and the force and moment equations 

(2. 58)— (2. 66). Doing such, one finds that the auxiliary condition, (2.107), becomes 

(7T 21 d2 +7r 22 D) ^ +( “xz D2 +1T ^ D + + <V^° 2 +Tr 9.7$ = °» ( 2 - 112 ) 


22 


23 24' 


25 


26 27' 


where 


’21 * R a^ - C Y. > - C n. 1 


’22 * -< C „ +R a C Y >• 


’23 » -< C „ + R a C Y >• 

P P 

’24 * S9 o (C n + ft a C Y > - C9 o [ftfe + ft a< 6 ’ m i»- 

V V 

ir or = i CO + i S0 , 

25 zz o xz o 


= C SO - C CO - R (CO C v - SO C ), 
26 no n o a'oY oY 
p r r p 


and 


"27 " + C n +R a C Y >• 
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Similarly, (2. 108) becomes 


<*28° 2 + ”29 D)i: +(,t 30 D2 + (‘yy^ + " 3 2 D + "33 )9 = °> 


where 


"28 " W C Z . - “> + C9 o C Z. 1 - < CO o C m. * S9 o C m . >• 
w u u w 

tt = b [C 0 (c_ - jx) - S 0 c„ ] + so c - co c , 

30 a 1 o' Z. o Z . 1 o m. o m. 

w u u w 


TT on = R (CO C„ +S 6 C„ ) - CO C - SO C 
29 a'oZ o Z ' om om 
u w u w 


it = R (C0 C_ -S0 C 7 )+S0 C - CO C , 
31 a'oZ oZ om om 

w u u w 


77 32 _ R a C Z " C m ’ 

q q 


and 


"33' C V C m - Vz > - S9 ol ftfe 4 R a (B - mg)], 

w w 

Also, in the same fashion, (2.98) becomes 

< C D 2 4C, D)J - (i^D 2 - C D + S8 0 C >? 4 (» 19 D 2 4 ^ 20 D 

V V p V 

+ c £ $= o, 

V 

where 

tt = i S0 +i CO and tt b C CO - C SO . 

19 xx o xz o 20 i o t o 

r p 

Further, the end condition, (2. 109), becomes 


where 


( a£-) = (^d 2 +Tr g D)y + (" 9 d +* 10 )4> +( 1t 11 d +Tr 12 $> 

a 


(C -p) c c 

v v £ 

^7 = * ’ ^8 = - ’ ^9 ~ - * 

T T T 

o o o 


(2. 113) 


(2. 114) 


(2. 115) 
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’io * r I(6 - m ‘ e)C8 o - s8 o c y i - s vS> - c vS> • 


"ll»^< C8 o C Y - S8 o C Y>- 
T r p 

o 


and 


Yy / aX \ 

"lZ*T - + W ' 

T a 

o 

Also, the end condition, (2. 110), becomes 

8z ! 2*2* ~ 

<W> =<n 13 D + ’ , 14 D)X * ln 15 D ■ t,r 16 D)Z+(, 'l7 D+ ’l8 )? - 
a 

where 


( 2 , 116 ) 


Tr io = — [ce ( ce cv + S0 c ) - se ( ce + so c ) , 

13 * 1 o' o Z. o Z. o v o X. o X. 

T u w u w 

o 


tt = — se (ce c + se c ) - ce ce c„ + se c_ ) , 

14 * 1 o' o X o X ' o' o Z o Z ,J 

T u w u w 

o 


tr = — [se (se c - ce c ) +ce (ce c„ - se c„ ) - n], 

15 * 1 o' o X. o X. o' o Z. o Z. 

T u w w u 

o 


"16 * r-tsvs^x - C8 o C X > +C8 o ,C8 o C Z - S 8 o c z»- 
T u w w u 

o 


’i7 e r (c8 ° c z - S8 o c x>- 

T o « « 


and 


tt 1q = -i-ce (se c v - ce c 7 ) - (|f) . 

18 * o' o X o Z ' '9s 

T w w a 

o 
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(2. 117) 


And finally, the end condition, (2. Ill), becomes 

9# T 2 2 

(_ 0s’ ) s (1T 1 D + 7T 2 D) * +(7T 3 D +1T 4 D ) 2; + ( Tr 5 D +1t 6 )^» 

a 

where 


and 


TT 


TT . 


TU 


=— [ce (c ce +c se ) +se ( ce c „ + se c„ )-n], 

* 1 o' X. o X. o' o' o Z. o Z. 

T u w u w 

o 


= — [ce (c: co +c v so ) + se (Co c „ + se )], 

* o' X o X o' o' o Z o Z ,J ’ 
T u w u w 

o 


b -^-[co (co c -so c v )+se (co c n - so c„ )], 

* 1 o' o X. o X.' o' o Z. o Z. 

T w u w u 

o 


= —ice (ce c v - so c v ) +so (ce c„ - se c„ )], 

,1 L o' o X o X 7 o' o Z o Z /J5 

T w u w u 

o 


= -^ce c v + se c_ ), 

* ' o X o Z 
T 0 q Q 


'T-t C V C(, o C X + S(, o C Z ) + B-m-gl + (g) . 
T w w a 

o 


2. 6 The Transformation of the End and Auxiliary Conditions to the Cable Coordinates 
Note that the end and auxiliary conditions, (2. 112)-^2. 117), are expressed 

A A A 

in terms of the x’ ,y' and z' coordinates of the cable, and the x,y, and z 
coordinates of the body' s mass center. Thus, in order to apply these conditions 
directly to the cable equations, (1. 28), (1. 31), and (1.34), one must transform them 
to the cable coordinates £,y and l. Consider now the following transformation 
equations (see Fig. 8): 

x = Cc?£(l,t) - Sff&(l, t) + R CibCO, (2. 118) 

a 

y = y+RSi pCO, (2.119) 

a 
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H 


and 

z = C5t,(l,t) + £K£(l,t) - R SO. (2. 120) 

3 . 

Using the small perturbation assumptions for the cable (1. 14)-*(1. 17), and the 
small perturbation assumptions for the body, (2. 19)-*(2. 22), one obtains the 
transformation equations for the first order problem. Further, nondimensionalizing 
by equations (1. 24), (2. 46), (2. 50) and (2. 51), one has 

Dx = - (~-) S?d£( 1 , t) - 2R a Se Q D?, 


2* 2L ~ 2r a 2~ 

D x = -(-r^SaD &(l,t) - 2R S0 D 6, 
b a o 


( 2 . 121 ) 


Dy = (^)Dy(l.t) +2R a C0 o D0, 

D 2 y = (4~)D 2 y(l, t) + 2R C6 d\ (2. 122) 

D EL O 

Dz = (^)CSDC(l,t) - 2$ C0 D?f, 


2* 2L ~ a ~ 2 ~ 

D z = (-j^CaD ?,(l,t) - 2R a C6> o D 6. 

Note now that, at the attachment point, one has the following relationships: 


(2. 123) 


<S> 


and 


(2. 124) 


,9z x ,9x. dt, j 9z 9^ * 

W = ( 9? 9^ 1>t)+( 9i ) 9^ M) * 


(2. 125) 


a a a 

Multiplying (2. 124) by (9z/9s) and (2. 125) by (9x/9s) and adding the two, one 

el a. 

obtains a relationship for the cable slope in the two coordinate systems. Using the 
small perturbation relations, (1. 14)-*(1. 17) and (2. 91), and nondimensionalizing 
by (1. 24), one obtains the relationship in the following form: 
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L 


H,. f . ,ax. ,3*' pi si' 

a a a a 


(2. 126) 


Now, when one uses equations (2. 121)— (2. 123) and (2. 126), the end and auxiliary 
conditions (2. 112)— (2. 117) may be expressed in the cable coordinates. First, 

(2. 121) and (2. 123) into (2. 113) give the auxiliary condition: 


where 


( * 45 D2 + * 46 D )£( M > + (* 47 d2 + Tr 48 D + * 33 )? = 0 , 


*45 ^ <f )( *30 CS - *28^’ 


*46 s <^><*31^ - *29^’ 


it.. = i - 2R (v S0 + tt„„C 0 ), 
47 yy a' 28 o 30 o' 


and 


(2. 127) 


it = v - 2R (ttS0 +ttCA ). 

48 32 a' 29 o 31 o' 

Similarly, (2. 122) into (2. 112) gives an auxiliary condition: 

( ”49 E2 + "50 D);?<1 '‘> +( -‘x Z d2 +r 23 D+ ’'24> ? +(,I 51 d2 + " 5 2 D 

+ " 27 $ = o, 


(2. 128) 


where 


,2L. 


,2L. 


*49 ~ * b**21’ *50 ~ ( b^*22* 


tt = 2R ce Tr +tt and Tr = ir +2R ir C 9 . 

51 a o 21 25 52 26 a 22 o 

Further, equations (2. 122) into (2. 114) give another auxiliary condition: 

<\3 d2 + ”54 D)5(1 '‘ ) + '-‘xx d2 +C i D - S9 o C i ft 

P v 


+ < ’65 D + *t 6 D + C 2^’° 

V 


(2. 129) 
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where 


*53 * <f> C *. • *54 ■ ‘iK • 

V V 


IT = ir +2R ce , and tt = tt +2R CB . 

55 19 a o 56 20 a o 

Now, an end condition is given by (2. 121), (2. 123), (2. 116), and (2. 117) into 

(2. 126). This becomes 

!§<l,t) - OT 40 D 2 + » 41 D)?(l,t) +(%2 D 2 4* 43 D 41 T"iS, 

where 

*40 * <*15^- *13^» - <f> ^S 05 ' 'I s ® 1 ’ 

a a 


(2. 130) 


2L. r .9x. 


,9z. 


*41 * ( T7 H W ('is 01 - *14 &) - W ( *4°* - *2 a,,) ' 


*42 * <’l S<> „ **3 CS o) - <t> <*13 S8 o + *15 C8 o>1’ 

a a 


",0 = [(f 2 -) - 2R (tt S0 +f, „C0 )"> - (~) /tt - 2R (ttS0 

43 lv 9s v. 17 a' 14 o 16 o'J v 9s X 5 a v 2 


9z. 


and 


*4 C<, <> 


,9x ,9z. 

*44 * W *18 ' <K» *6' 
a a 

Also, an end condition is given by (2. 122) into (2. 115): 


(V 

ga.f) - <* 57 D 2 4- 68 D)y 4 ,„ 9 D 4„ 10 )? 4( %9 D 2 4. 60 D 4H 12 )?, ,2. 131) 


57 


where 


,2L V .2L. 

*57 5 <T>*7- *58* < “b ) *8- 
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TT 5 ; 2R ce TT and tt = tt +2R C0 tr 
59 a o 7 60 11 a o 8 


Finally, to complete the set of end conditions, assume that the cable is fixed at the 
origin of the x,y, i coordinate system. This gives 


U 0,t) = 0, (2.132), y(0,t) = 0, (2.133), and C(0,t) = 0. (2.134) 

Thus, (2. 127)-*(2. 134) give the end and auxiliary conditions for the cable equations 
(1. 28), (1. 31), and (1. 34). The nature of the solution is such that initial conditions 
neqd not be specified, that is, the general solution is sought. Therefore, the 
problem statement is now complete. 
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3. THE SOLUTION OF THE CABLE-BODY EQUATIONS 


3. 1 A Discussion of the Equations 

Note the very important fact that the equations uncouple into two separate 
problems. The cable equation, (1.34), along with the end conditions, (2. 130) and 
(2. 134), and auxiliary condition, (2. 127), constitutes a complete problem for the 
general solution of £(s,t) and 0(t). Similarly, the cable equation, (1.31), with 
the end conditions, (2. 131) and (2. 133), along with the auxiliary conditions, (2. 128) 
and (2. 129), gives a complete problem for the general solution of y(s, t) , $>(t) and 
% (£). Physically, this means that the first order problem uncouples into two 
distinct modes: lateral and longitudinal motions. Such uncoupling is, in fact, 
observed by experiment (Chapter 4 and Ref. (1)). Further, these motions have 
certain analogies to uncoupled airplane motion, although the comparison must not 
be carried too far, since the two mechanical systems are fundamentally different. 
Each of these problems is now treated separately. 

Note finally that, consistent with condition (1. 15), the £(s,t) coordinate is 
of no significance in the first order problem. 


3,2 The Longitudinal Solution 

The longitudinal problem is described by the following equations: 
Cable Equation: 


D2? - 62 5 +k 3 D?+k ^i =0 - 

8s 


End Conditions: 


(1.34) 


= (* 40 D 2 +ir 41 D)£(l,t) + (u 42 D 2 +tt 43 D +tt 44 )2, (2. 130) 

and 

£(0,C) = 0. (2.134) 

Auxiliary Condition: 

(1T 45 D2 +7r 46 D ^ (1 ’ f) +(TT 47 d2 +Tr 48° +1T 3 3^ = °* < 2 ' 127) 
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These equations are linear, and have constant coefficients. Therefore, in the 
spirit of linear differential equation theory, the solution is assumed to be of the 
following form: 

A A 

£(s,t) = Z(s)e°^ and 0 = ©e°*, (3. 1) 

where © is a constant. Now, substituting this into the partial differential equa- 
tion, (1.34), one obtains an ordinary differential equation for Z(s): 


. 2 , , dZ -2d Z 

(U + kg(7)Z + k 4 ^-C — =0. 

ds 


(3.2) 


Also, by ordinary differential equation theory, the general solution of this is 

Z(s) = Z e (A+ 0)8 + Z e (A " fi) ®’ (3 ‘ 3) 


where 


A == — r and 
A 2 
2C 


»■[ 


k 4 (ff 2 +k,u )- 1/2 

*4 A 2 

4C C 


: 1 ■ 


(3.4) 


Now, by the end condition, (2. 134): 

Z(0) = 0=Z +Z . 

°1 °2 

Thus, (3. 3) becomes 

„ As. ns -Os. 

Z(s) = Z Q e (e - e ). 

Using this, and substituting (3. 1) into the end condition (2, 130), one obtains 

e V - e _n ) *V A ' Z o + <V’* 

+ ir 43 ff + ir 44 > 0 = °- 

Further, (3.5) and (3. 1) into the auxiliary condition, (2. 127), give 
eV - e‘V 4R a 2 + * 46 a)Z o + W +1T 33 )e = °* 


(3.5) 


(3.6) 


(3.7) 


45 46 ' 0 

Equations (3. 6) and (3. 7) are two linear homogeneous equations in Z q and 0. 
Thus, it follows that an equation for a (characteristic equation) may be obtained 
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by putting these two equations into a determinant, and setting it equal to zero. 
Doing this , one obtains 


<*,,/ + ' r 4i t ' - A - n c°th n), (tr 42 cr 2 + » 43 <r 




( V + V 


where 


= 0, 


(3. 8) 


coth 0 = +e 

(e n - e"°) 

Note that a is, in general, complex; that is, 

• , ,.l/2 

(T=cr r +)(7y where 3 = (-1) 


(3.9) 


Thus, the characteristic equation, (3. 8), is a complex transcendental function in 
a complex variable. To facilitate finding the roots of this, one expands it into 
two real characteristic equations in two real variables, cr^ and 17.. To this end, 
note first the expansion of fj coth n : 


n = n. + jo., 


where 


h 2 h 2 
n = — cv, n. - -%v, 

r 1 a 

c 1 c 


k 4 2 2 2 2 l/k 
h 2 s [( ^2 + ^r " a j + W + < 2tr r ff j + Vj* J ’ 


and 


(3. 10) 


—1 k 

V , JS |-t<2c V r J + Vj >/( ^| +< r - VVi 11 ' 


Thus, 


ncoth n= H 3 +jH 4 , 


(3. 11) 
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where 


2 2 
ft coth 0 (1 +cot ft.) - ft. cot ft.(l + coth ft ) 
u = r r' 3 3 3 r 

3 " .. 2 2 

coth ft + cot ft. 
r 3 


and 


H 4 2 * * S 


ft. coth ft (1 +cot ft.) + ft cot 0.(1 - coth 0 ) 
J 3 £ i r_ 

2 2 
coth 0 +cot ft. 
r 3 


Now, this and (3. 9) are substituted into the characteristic equation (3. 8). 

Expanding this, and separating into the real and imaginary parts, one obtains 

two simultaneous real characteristic equations in two real variables, and 

<7.. These are 
3 


and 


where 


r' r’ j' 


3 r 3 


E 1 * V^r - "?> Vr ' (A + H 3 )> 

E 2 e ^oVj +lt 41 ff j * *V 
2 2 

E„ s u fcr -o\)+tr a + ir_,_,. 

3 42' r j ; 43 r 44* 

E . = 2ir <7 a. +ir or., 

4 42 r j 43 3 

2 2 

E = tt (a - a. ) 4ir a , 

5 45 r j 46 r 

E„ s 2 tt cr (7. + 1 T o., 

6 45 r 3 46 3’ 

2 2 

E = tt (a - cr ) + tt a + n- # 

7 47 v r 3' 48 r 33’ 


+ E 6 E 4 = 0 

(3. 12) 

- E 6 E 3 - «• 

(3. 13) 


and 


E d ==■ 2tt a cr. 4 tt < 7.. 
8 47 r 3 48 3 
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An electronic computer is used to solve equations (3. 12) and (3. 13). This is 
explained in detail in Section 3. 4. 


3. 3 The Lateral Solution 

The lateral problem is described by the following equations: 
Cable Equation: 

.2- 


8s 


End Conditions: 


#N/ A 


and 

y(0,t) = 0. 
Auxiliary Conditions: 


and 


V 2 + ’ , 50 D tf (M ‘ ) ^-‘xzH 2 +ir 23 D+,r 24 )? +(,, 51 d2 + ’' 5 2 I> 
+ W 27^ = ° 


<’53 d2 *” 54 D >y<M> *c t D - sec^ + <-> 5 ,V 

*\6 D * C l 

V 


(1.31) 


g(i,t) = (u 57 d 2 +1T 58 D )y +( 1t 9 d+1t 10 )^ +(^ 59 d 2 +-r 60 D + " 12 $ (2. 


131) 


(2. 133) 


(2. 128) 


(2. 129) 


As in the longitudinal case, these equations are linear, and have constant coef- 
ficients. Therefore, the solution is assumed to be of the form: 

A A *. 

y(s,t) = Y(s)e Xt , $ = tfe**, and $ = (3. 14) 

where ^ and 3> are constants. Now, substituting this into the partial dif- 
ferential equation, (1.31), one obtains an ordinary differential equation for Y(*5): 


< x2+ V> Y+ sf- £2 fi=o- 

ds 

Further, the general solution of this equation is 


(3. 15) 
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where 


Y(s) = Y e (r+A) *+Y e (r - A)i , 
°1 °2 


k ? ( x2+k 6 A ) 


and A = — - + 


“ *2 -2 -2 * ' ' 
2C )_4C C J 

The end condition, (2. 133), gives 

Y(0) = 0 = Y + Y . 

°1 °2 

Thus (3. 16) becomes 

Y(s) = Y Q e rS (e A ® - e -A ®). (3. 

Using this, and substituting (3. 14) into the end condition, (2. 131), one obtains 

e r (e A - e‘ V 5? X 2 + » 5( X - r - 

(e -e n 


* < V + "lO*® + (,r 59 X + "eO* + ,t 12 ) * " °- 


Further, (3. 18) and (3. 14) into the auxiliary conditions, (2. 128) and (2. 129), give 


S r (e A - e‘V 4 „X 2 + ’ r 60 x ) Y „ + <-‘ x2 x2 +r 23 X + V® 


+ (iT X +1T X+TT Yjr = 0 
1 51 52 27' 


(3. 20) 


e r (e A -e'V r A 2 +u r 1)Y +(-i X 2 +C X - SO C. )$ 

' 54 ' ft v vv I n # . 


53 54 o xx 


4 <’55 X +1T 56 X +< “i W= °- 

V 


(3. 21) 


Equations (3. 19), (3.20), and (3.21) constitute a set of three linear homogeneous 
equations in Y q , % and <£. Thus, as in the longitudinal case, a characteristic 
equation for X is obtained by putting these into a determinant, and setting it 
equal to zero. One then obtains 
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( "57 X2 - r ' A COth A >- <V +, 'lO )l + ”60 X + ”l2 ) 

( " 49* 2 +,, 50 X) - + ” 23 X+ ” 24 ) ’ + " 52 X+ ’' 27 > =» ( 3 . 22 ) 

+ ”64 X >' <-‘xx X2 +C jt X - Se 0 C i »■ + ” 6 6 X +C 1 > 

P V V 


where 


coth A s 


(e A 4e -A ) 

(e i -e- A ,‘ 


Note that A is, in general, complex; that is 


X = X^ + j?y where j = (-1) 


(3. 23) 


Thus, the characteristic equation, (3.22), is a complex transcendental function 
in a complex variable, X. As in Section 3. 2, the equation is expanded into two 
real characteristic equations in two real variables, X^ and X^. Note now the 
expansion of A coth A: 

A=A r +jA., (3.24) 


h l h l 

A = — Cv, A. = — Sv 
r » i » 

C 1 C 


k 7 2 2 2 2 lA 

h, = [( — ^r + X^ - X; +k X ) +(2X X. + k X.) ] , 

1 4C 2 r * 6 r r J 6 j 1 


-1 k 

u l2£_ [(2 X X. +k X.)/(-l- + X 2 - x 2 +k X.)]. 
2 lv r j 6 V 4 g2 r J 6 j n 


A coth A = Hj + jH 2 , 


(3. 25) 


where 


J 


H l 3 


A coth A (1 +cot^A.) - A. cot A.(l - coth^A ) 

r* v v i i t' 


J1 l 

2 2 
coth A +cot A. 
r 3 


_L 


and 


A.coth A (1 +cot^A.) + A cot A.(l - coth^A ) 
it _ 3 r v 3 r 3' 

n 2 ~ 2 2 

coth A +cot A. 
r 3 

This and (3. 23) are substituted into the characteristic equation (3. 22). Expanding 
this, and separating into the real and imaginary parts, one obtains two simul- 
taneous real characteristic equations in two real variables, X^ and X^. These 
are 

F r (X r’ V = A 1 A 19 ~ A 2 A 20 - A 11 A 21 + A 12 A 22 + A 17 A 23 ' A 18 A 24 = ° (3 


and 


where 


F J^r*V A 2 A 19 +A 1 A 20 " A 12 A 21 “ A 11 A 22 +A 18 A 23 +A 17 A 24 °’ (3 ' 


A l S ’'57< X ?-^ + ’'58 X r- r - H l- 

A 2 ' 2 ,r 57 X r X j + "58 X j * H 2’ 

2 2 

A„ = -i (X - X. ) + ir X +r „ 

3 xz v r j' 23 r 24* 

A =-2i XX+ir X, 

4’ xzrj 23 j* 

A = IT ,(X 2 - x 2 ) + it X + tt , 

5 51' r j 7 52 r 27’ 

A. = 2tt X X. + tt X 

6 51 r j 52 3 

A =-i (X 2 - X 2 ) + C X -SO C. , 

7 xx v r j l r o l 

P v 

A 0 = -2i X X. + C X. 

8 xx r 3 jt p 3 > 

A 9 sTr 55 (X ? _X ? )+Tr 56 X r +C X * 

J v 


A „ = 2 tt XX+tt X, 
10 55 r j 56 j’ 


.26) 

27) 
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and 


A = tt (X 2 - X 2 ) +tt X 

11 49' r j' 50 r 

A,„ = 2tt X X. +u X 

12 49 r j 50 j’ 

A 13- Vr + ”l0- 

A X4 * Vj> 

2 2 

A „ _ = ir „(X - X ) +tt X +tt „ . 

15 59' r j' 60 j 12’ 

A„„ = 2tt X X. +tt X . 

16 59 r j 60 j’ 

* .-.2 ..2 

A„„ = tt (X - X,) + tt X , 

17 53' r y 54 r’ 

A 18 S **53 V] +U 54 X j’ 

A 19 S A 3 A 9 + A 6 A 8 “ A 4 A 10 " A 5 A 7’ 

A 20 S A 3 A 10 +A 4 A 9 ” A 5 A 8 " A 6 A 7* 

A 21 s A 13 A 9 +A 16 A 8 " A 14 A 10 " A 15 A 7’ 
A 22 S A 13 A 10 +A 14 A 9 “ A 15 A 8 " A 16 A 7’ 
A 23 55 A 13 A 5 +A 4 A 16 ~ A 14 A 6 ~ A 3 A 15’ 


A 24 - A 13 A 6 +A 5 A 14 “ A 15 A 4 A 16 A 3’ 


An electronic computer is used to solve equations (3. 26) and (3. 27). This is 
explained in Section 3. 4. 


3. 4 The Computer Solution of the Characteristic Equations 

Since equation set (3. 12) and (3. 13) and equation set (3. 26) and (3. 27) are 

mathematically similar, that is, both sets are two simultaneous nonlinear 

transcendental equations in two unknowns, the method of root extraction applies 

to both cases. This method is a roots locus plot, such as used by Dugundji and 

Gareeb (Ref. 3). For example, consider solving for the X roots. X^ and X. 

are systematically sequenced through a range of values. For each of these 

values, F and F. are calculated. Now, for each X ,X, pair for which either 
r j r ] 
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F or F. equals zero, this X , X. location is marked on a X , X. coordinate 
r j r j r j 

system (Fig. 9). Thus, after sequencing X and X, through their full range 

^ 3 

of values, one obtains a series of these zero points through which one may 

draw F = 0 and F. = 0 curves (Fig. 9). The intersection of a F =0 curve 
r j ' r 

and a F^ = 0 curve defines a X root on the coordinate system. 

Now, an electronic computer was used to find F^ and F. for sequenced 
values of X^ and X^ (similarly G r and for c and cr^) ; but the plotting 
and root extraction was done by hand. This was done in order to keep the 
computer run time to less than one minute. However, the hand plotting was by 
no means difficult and gave roots in less than five minutes. Also, the plots had 
the virtue of showing trends, when compared with other plots in a series of runs. 
Examples of roots locus plots are shown in Figures 18 and 19. 
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4. A COMPARISON OF THEORY WITH EXPERIMENT 


4. 1 The Test System 

For the test system, akite, as shown in Fig. 21, was chosen. This 
choice was based on the fact that a kite contains all of the essential features of 
the theory' s physical model and is much more general than a nonlifting design, 
such as tested by Etkin and Mackworth (Ref. 5). Further, the kite was inex- 
pensive to construct, easy to test, and perfectly suitable for the Stanford low 
speed open throat wind tunnel, which has a test section diameter of 7 ft. By 
virtue of the tunnel design and low wind velocities, it was possible to study the 
system's unstable motion without great risk of damaging it. 

The system, however, has certain limitations. Aside from being repre- 
sentative of only one cable-body application, wind tunnel test section size limited 
the cable length such that cable dynamic terms were very small compared with 
body terms. To offset this, a certain amount of outdoor testing at long cable 
lengths was attempted. Although outdoor experimentation is difficult because 
of atmospheric vagaries, as discovered by Bryant and coauthors (Ref. 1) and 
NA SA researcher Tracy Redd (Ref. 15) with his towed balloon experiments (fig. 17) > 
reasonable information was obtained, and the sum total gave a good spectrum of 
the system's stability characteristics. 

As for the body itself, note from Fig. 20 that its design was kept purposely 
simple to facilitate calculation and measurement of its properties. By virtue 
of this , the moments of inertia were evaluated by a combination of theory and 
compound pendulum tests. Their values with respect to the reference axes 
(Fig. 20) are 

(I ) = 2. 57 X 10~ 2 slug-ft. 2 , (I ) = 2. 99 X lo“ 2 slug-ft. 2 , 

xx'o yy 0 

(I ) = -. 177 X 10 -2 slug-ft. 2 , and (I ) = 5. 50 X 10 -2 slug-ft. 2 . 
xz zz 

o o 

Further, the stability derivatives were calculated according to Campbell and 

McKinney (Ref. 2), Etkin (Ref. 4), and Purser and Campbell (Ref. 14). A 
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certain amount of difficulty was encountered in calculating the contribution of the 

"V" tail to the stability derivatives. The theory of Ref. 14 is primarily for tails 

with dihedral angles of 30 degrees or less. Moreover, no information on a "V" 

tail' s contribution to Cy and C« was available. However, for lack of a 

P *P 

better method, (C Y ) ,(C« ) and (Cy) were found from References 2 and 
p T *PT « T 

14; and (Cy ) and (C* ) were calculated based on assuming the section lift 
P t P T 

curve slope to vary elliptically over the tail' s span. An estimation of the error in 

(C Yg) » (0^)^, and (Cy^^ is given in Section 4. 4. The primary motivation for 

using a "V" tail, as opposed to a more conventional design (for which there is 

much more information), is that a "V" tail is much less susceptible to damage 

in the case of the body's tumbling Further, a tethered lifting body often flies fully 

stalled in certain conditions, for instance, during launching, and experience shows 

that a "V" tail gives superior directional stability for this situation. 

In order to facilitate studying the tail' s contribution to the lateral stability 

derivatives, in particular, the effects of (Cy.) ,(C«„) , and (Cy ) , the 

x p t ^ A P T 

author wrote the lateral stability derivatives as functions of these terms. Note 
also that these stability derivatives are with respect to the wind reference axes 
of standard aircraft practice (see the Appendix). 

The values calculated for the stability derivatives are 
<c x ) =0, (C z ) =0. (C m ) -o, (C x ) =0, 

U u U U 

O 0 O O 

< C Z.> = 0 ' < C m.» '°-< C X> --223C l . 
u u a o 


(C„ ) = -4. 12 - C n , (C ) = 2. 61, (C ) = 0, 


a 


m 


a 


X. 

O' 


(Cry ) =-1.58, (C_ ) =-3.41, (C x )=0 


Z. 

a 


m. 

a 


(C ) =5.42, (C) =-17.10, (C ) =-.75, 


m 


P 


(C v ) = -.05, C = (C ) +(C ) , 

6 (3 6 0 

p w H H T w 
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(C ) = 45(C ) + . 225(C ) + . 02CT , 

B B y B w 

y o H T p w 

(C ) = -. 154 +. 336a , (C ) = -.127 -.07C_ , 

P T *0. 

<C. ) = (C ) + (C ) , (C ) = 2(C ) , 

P P „ P r n p 

o k T p w o k o 

(C ) = 90(C ) + . 45(C ) - . 023 C 2 - . 30(C. ), 

n n 0 Li u 

r P rr> P W W 

o T w 

(C ) = 90(C ) + .45(C ) +. 26C +.127, 

L Xry Li 

r S 6 w 

o r T ^ w 

(C y ) = -.207 +.907a b , (C y ) =-.225, 

P T P W 

(C y ) = (C y ) +(C y ) , 

P o p T P w 

<C n ) = -. 45(C y ) + . 23(C y ) -.06 C l +9.0(C ), 

P Pm P W W 

o T w a 

(C^ ) =-.379 -.13C d -. 0012(3. 78 - 16. to b ) 2 , 

p w 

o 

(C ) = 0, (C ) = 0, and (C ) =0, 

P P P 

r o y o r o 

where = the body' s fuselage angle in radians (see Fig. 10). For these 

equations, the subscript ( ) refers to the wing. Also, the cable used was 

w 

stranded nylon with the following properties: 

p = 1. 78 x 10~ 5 slug/ft. , R = . 0023 ft. , 

and according to Hoerner (Ref. 7), 

C = . 035 and K = 1. 15. 
a 

o 


Now, the equilibrium configuration of the system for given wind speeds, 

A 

U, and cable attachment point, R , Is specified by the quantities a, T , and 

& o 

Values for these were obtained directly from the experimental investigation 
of a large number of equilibrium situations. A beam balance was used to 
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measure the vertical component of the cable tension, T Q Sin a, and a 3 degree 
inclined manometer measured the velocity, U. Also, the angles a and 0 
were found directly from photographs of each test (Fig. 21). From these data, 

a 

values for T were calculated from 
o 


T = 
o 


2(T Bin a) 
o 

t T 2„ . ~ 

pU S sin a 


(4. 1) 


Graphs of T ,a, and 0 vs. U are given in Figures 22, 23, and 24, and the 
o o ' 

estimated error in these values is given in Section 4. 4 and shown in the graphs. 

Also required, for the stability derivative equations, are the equilibrium 
values: 0^,(0^), and These are found by a 

combination of theory and the experimental equilibrium quantities. Cj^ is 
given by 


C =T o sS-mg, 
o 


(4. 2) 


and Cjj is calculated from 
o 


C D 'V*' 
o 


(4.3) 


was measured directly from the photographs; and, from the geometry of the 
body, one has that 


a - of, +4. 5 . 
w b 


(4.4) 


Now, in order to obtain the wing properties, the author used the experimental 
curves by Pinkerton and Greenberg (Ref. 13) for a rectangular wing with the 
same airfoil, but of aspect ratio 6. The angle of attack of a wing of aspect 

ratio 6 such as to give the same C as a wing of aspect ratio 4. 83 is given by 

JL 


(<1+(C L ) /6TT) 


a 


00 


01 6 (1+(C T ) /4.83irf4. 83 ' 947q V 


(4.5) 


a 


oo 


So, this equation and the experimental curve of versus a gives Cj^ for 

each a value. Similarly, the angle of attack for a wing of aspect ratio infinity 
w 
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(4.6) 


is related to a by 
w 


ot = . 727a . 
00 w 


This equation, with the experimental curve of C<j w versus gives Cj w 

for each a w (and thereby Cj^) value. Further, the slope of the curve at that 

point gives (C^ ) . Finally, C D is calculated from 
w a w 


C,. = C , + . 066C 2 . 

Da L 


(4.7) 


w w w 

It is important to note that for all equilibrium quantities except a, the 
effects of the cable weight and drag are ignored. This is by virtue of the fact 
that a cable length of 4 feet or less was used during the equilibrium measure- 
ments, and cable weight and drag were very much less than cable tension and body 

drag for the wind velocities, U, considered. This gave an essentially 

straight cable which directly yielded 


,d\ , ,8z, „ 

(— ) = Ca and (— ) = Sa (4. 8) 

a a 

Also, since most stability measurements were made at a cable length of 4 feet 
or less, the assumption of constant equilibrium parameters for a given U and 
R , and varying L, is very good. However, for the case where L = 100 ft. , 

Si 

note again that the cable weight is still much less than for the range of U 

and R values considered, that is, 
a 


pgL = 0[10 _2 ]lb. and T q = 0[l]lb. (4. 9) 

But, the cable drag is of the same order as the body drag: 

pKRU 2 L = 0[10 -2 ] and D = 0[10~ 2 ]. (4. 10) 

Its only contribution, though, is to give the cable a bow, and to lower the cable 
mean angle, a. The cable slope at the attachment point remains unchanged. 
Thus, the equilibrium quantities for the L = 100 ft. case are the same as for 
L = 4 ft. or less, with the exception that a is different. 

In order to find the change in a, note first that (4. 9) and (4. 10) satisfy 
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conditions (1. 16). So, the cable bow is shallow. Also, for this particular model, 
the a value is very high (70° or more), so that one has for the angle decrease, 
Ac?, 



o 


This now completes the set of parameters and equations necessary to describe the 
equilibrium configuration of the system. These parameters, along with the 
inertial properties, stability derivatives, and characteristic dimensions give all 
the information needed by the theory to define the stability roots. 

4. 2 The Stability Tests 

For a given free stream velocity, U, and attachment point position, 

R , the stability quantities measured were the cable' s critical length and the 

2r 

system' s lateral and longitudinal oscillations. The cable' s critical length is 
that value of L at which the system went unstable. This appears to be an 
important phenomenon of cable-body systems, and has been observed by Bryant, 
et al. (Ref. 1) and Etkin and Mackworth (Ref. 5). As for the oscillations, the 
properties recorded were the frequency and qualitative damping of the lateral 
and longitudinal motions. Finally, a very qualitative measure of the system's 
stability at long cable lengths was made. 

The methods of testing were very straightforward and direct, which was 
one of the virtues of using a low speed kite system. Critical cable length was 
measured in the wind tunnel by slowly unreeling the cable until unstable oscil- 
lations occurred. The cable was then marked, and the critical length, L , 

cr 

was directly measured. To produce the lateral oscillations, '’'■jgjjp* and the 
longitudinal oscillations, °j e » at L < L cr > the cable was perturbed a given 
amount by hand and then suddenly released (Fig. 12). For this system, it was 
very easy to produce almost pure lateral or longitudinal motions by this method. 
Frequency was then measured by using a stopwatch and counting integral numbers 
of periods. Damping measurement was somewhat more qualitative, and involved 
counting the number of periods until equilibrium was essentially reached. 
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Stability studies for the long cable case involved outdoor flying, which thereby 
restricted quantitative measurements. However, it could be determined 
whether the system were indeed stable or not, and if so, 

it was. Moreover, the system was stable during rapid unreeling of the cable. 

So it was possible to test unstable motions for values of L much larger than L 

Now, for those properties for which quantitative measurement was 

possible, L cr> ^j ex p» and cr j ex p» the repeatability was very good, and the 

values given (Figs. 25 and 26) are based on the average of several runs. The 

largest estimated error in L is 1 . 2 ft. ; and for X* and crs , the 
& cr Jexp Jexp 

largest estimated error is ± . 01. Further, the maximum estimated error in U, 
for the velocity range considered, was 6%. Somewhat large tunnel turbulence 
precluded an accurate damping measurement. This was the reason why a movie 
camera was not used. It was actually more meaningful to measure the damping 
as described, and to make a qualitative judgment. Insofar as comparisons of 
damping from one test to the other were concerned, this was quite adequate. 


cr 


4. 3 The Computer Examples 

For the experimental cases, the wing lift coefficient, C is no less 
than .70. Moreover, for a wing operating at a high Cj^, the stability deriva- 
tives may be strongly modified by partial flow separation. Thus, since the calcu- 
lation of the wing' s contribution to the stability derivatives is based largely on 
the assumption of attached flow (Ref. 2), the lowest value of U that was 
selected for the computer examples is 25. 0 ft. /sec. Nevertheless, was 

still higher than desired, and it is felt that this introduced a certain amount of 
error in relating the theoretical results to the experimental results. This is 
discussed further in the next section. 

For the lateral stability study, the particular examples chosen are 
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R = . 948 ft. 
a 


U = 25. 0, 28. 0, 30. 0 ft. /sec. / R = 1. 09 ft. 


R = 1. 24 ft. 
a 


L = 1. 00 ft. 
L = 2. 00 ft. 
L = 3. 00 ft. 
L = 1. 00 ft. 
L = 2. 00 ft. 
L = 3. 00 ft. 
L = 1. 00 ft. 
L = 2. 00 ft. 
L = 3. 00 ft. 


U = 30. 0 ft. /sec. , R = 1. 09 ft. , L = 100. 00 ft. 

cL 


For the longitudinal stability study, the particular examples chosen are 

r R = 1. 09 ft. , L = 2. 35 ft. 

U = 25. 0, 28. 0 ft. /sec. / a 

1 R = 1. 24 ft. , L = 1. 76 ft. 

The equilibrium data for these examples is obtained from Figures 22,23,24, 

and 27, and equations (4. 2)-*(4. 7). The values obtained are tabulated below: 

U | 25. 0 I 28. 0 | 30. 0 

(ft. / sec. ) 


25. 0 

28. 0 


. 948 

1. 09 

1.24 

. 948 

1. 09 

1.24 

. 948 

.972 

.650 

.460 

.938 

.592 

.407 

.903 

-. 541 

-.532 

-.497 

-. 576 

-.562 

-.523 

-. 592 

. 134 

. 040 

0 

. 100 

.020 

-. 028 

. 078 

1.21 

.898 

.700 

1. 13 

.790 

.600 

1. 07 

. 176 

.111 

.091 

. 151 

.099 

. 086 

. 138 

1.42 

1. 03 

.94 

1.31 

.97 

.82 

1.27 

. 159 

.088 

.075 

. 137 

.079 

. 059 

.110 

. 026 

. 018 

. 016 

. 024 

.017 

.014 

. 023 

. 0022 

. 0019 

.0017 

..0021 

. 0018 

. 0015 

. 0021 
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For L less than 5 ft. , 


u 

(ft. /sec. ) 

25. 0 

28.0 

30.0 

R a (ft. ) 

.948 

1.09 

1.24 

.948 

1.09 

1.24 

.948 

1.09 

1.24 

a(rad) 



1.364 

1.412 

1.408 

1.358 

1.422 

1.406 

^9 


and for L = 100. 00 ft. , the previous a values and equation (4. 11) give 


U 

(ft. /sec. ) 

30.0 

R a (ft. ) 

1. 09 

a (rad) 

.958 


4. 4 Computer Results and Comparison of Theory with Experiment 

Using the method described in Section 3.4, one directly obtains values 

of and a r ,cr j ^ or com P u ter examples listed in Section 4.3. Further, 

for the lateral case, the range of L values went through a critical length, 

L , at which the system had theoretical neutral stability. In order to obtain 
cr 

L cr and Xj f the following technique was used. On a X^, X. versus L 

coordinate system (Fig. 13), points of X^ and X^ are plotted and curves are 

drawn through them, thus giving a X^(L) curve and a X (L) curve. Further, 

the intersection point of the X^(L) curve with the L axis gives L^. Finally, 

the X. value corresponding to L (X.) is X^ 

J cr' j' J C r 

Now, for the lateral case, the theoretical and experimental results are 
listed on the next page. 
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And, for the longitudinal case, 


u 

(ft. /sec. ) 

R 

a 

(ft.) 

L 

(ft.) 

£7. 

J 

cr. 

J exp 

a 

r 

qualitative 

experimental 

damping 

25. 0 

1.24 

1.76 

.472 

.437 

-. 072 

lightly 

1. 09 

2.35 

.465 

.458 

-. 110 

moderately 

28. 0 

1. 24 

1.76 

.450 

.435 

-. 058' 

lightly 

1. 09 

2.35 

.447 

.442 

-. 095 

moderately 


Notice that the theoretical results are somewhat different from the 


experimental values. L is uniformly less than L cr and Ai is con- 

ci* sxp Jcr 

sistently greater than (A; ) . in order to investigate the reasons for this, 

•' cr exp 

the author studied the effect of the input data on A^ and A^. Considering first 
the equilibrium values, the author found the significant parameters to be T q , 

6 , U, and a. So, starting from the U = 30. 0, L = 2. 0, R =1. 09 case (as a 

O cl 

representative example), a finite difference study using the computer gave 
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ax 

ax 

ax 

r 


• 25 ' a# 

8T 

80 ~ 
o 

o 



ax. 

dx . 

ax. 

L = 

. 10, — ^ = . 

19. i = 

at 

* dO 

0 

’ au 


ax 

01 , ~£~~ = • 06 , 
ax. 

I, and . 03, 

Also, the estimated error in the equilibrium values gives 


AT =.08, A0 =.06, AU =1.8, and Av = . 06. 


So, from 


AX = 
r 


dX 

r * 


dX 

r . ^ 


ax 

r „ 


ax 

r ~ 

AT 

* o 

dT 

o 

+ 

- — AO 

80 O 

0 


au A 




and 


AX. = 

J 


one obtains 


ax. 


ax. 


ax. 


ax. 

— ^-AT 

at ° 

O 

+ 

aa 1 Ad o 

O 

+ 

w au 

+ 



AX = . 034 and AX. = . 021. 
r J 

Thus, whereas errors in the equilibrium values have a small effect on X , 

J 

they have a significant effect on X^, and thereby, L . 

A second source for error was the stability derivatives. One problem, 
as mentioned before, was that the wing was operating near stall for the test cases 
considered. This was unanticipated as it was originally planned to test at a much 
larger U, and, consequently, at a much smaller Cp^. But the body' s structure 
was dangerously strained at U greater than 31. 7 ft. /sec. , so the test speeds were 
lower. Further, the airfoil chosen has very gentle stall characteristics; thus the 
body could fly well, even though a significant part of the wing' s flow was separated. 
Under such conditions it is not only difficult to assess the wing' s contribution to 
the stability derivatives, but one also encounters the question of nonlinear effects. 
Nevertheless, for lack of other methods, the wing's stability derivatives were 
calculated using Ref. 2. Note, however, that one of the wing' s stability 
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derivatives that is most affected by partial stall is (C n ) . For the U = 30. 0, 

P w 

R = 1. 09, L = 2. 0 case, 
a 


9A 

r 

9(C ) 

n 

P 


9\ 


= 3, and 


J 


w 


8<C n» 

P 


= - 2 . 0 . 


w 


Thus, even though it is difficult to give a value to the error in (C n ) , it is seen 

P w 

that relatively small changes in its nominal estimated value of . 03 can make 
large changes in A.. 

Another source of error in the stability derivatives was the "V" tail. 

Its area is large (S y /S = . 48), and thus it gives a major contribution to the 

stability derivatives. Unfortunately, the study of wings of very large dihedral 

has never been a very popular research topic, and the only report available for 

estimating some "V" tail properties is Ref. 14. However, it was desired to 

use a "V" tail for this body, for the reasons mentioned earlier; thus the tail' s 

stability derivatives were calculated as best as was possible. The significant 

tail stability derivatives are ( C Yp^ T » ^p^’ and ^ Cy p^x* For the * ong tail 

moment arm of the bodjr, the other tail stability derivatives, such as (C n ) , 

n P T 

come directly from the previous ones. Now, the estimated error in these 
values is 

A(C ) = .3, A(C ) =.04, and A(C V ) =.02. 

p T p T p T 

Further, the partial derivatives for the U = 30. 0, R =1. 09, L = 2. 0 case are 

Or 


9A 


9A 


9A 


9(C ) 

P. 


= . 10 , 


9(C ) 


= -. 20 , 


9A. 


9(C y ) 

P, 


9(C y ) 


= 0, 


9A 


= -. 16, 




T 

9A 


9(C. ) 

V 


= . 30, and 


9(C y ) 


= . 12 . 


So, from 
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ax 


dX 


ax 

d(C ) a(C y> 

S p T 

P T 

.1 

9 (C ) A(C £ > 
£ P t p t 


a<c y r ) * c v > 

P T T 

*r 

4- 


and 


M i = he \ a(C y > 

J 0(C Y„ ) Y p ^ 
P x K T 


+ ' *<%> + 

X ft Prr * P - 


one obtains 


A\ = . 042 and A*.. = . 062. 
r 3 

Thus, errors in the estimation of the tail stability derivatives are seen to have 

a significant effect on \ and L cr (through X f ). 

Note, though, that even with all of the sources of error combined in the 

worse possible way, the theoretical results should still be at least within the 

same order of magnitude as the experimental data. Indeed, this is not only so, 

but even better, for the error in Xj cr from experiment is no greater than 40 %, 

and for L , the error is no greater than 80%. More important, the theory 

predicts the essential features shown by the experiments. For instance, not only 

does it predict an L for lateral motion, but it shows the correct variation of 
cr 

L with U and R . Moreover, for the U range considered, both theory and 

Cl* a 

experiment show a decrease in Aj with U. Note, however, that the variation 

of (A* ) with R is difficult to assess within the error limits of the experi- 

Jcr exp a 

mental data; but the theoretical results likewise show no conclusive trend. 


Similarly, for the longitudinal motion, both theory and experiment show no definite 

variation of cr. with U and R , for the cases considered. But otherwise, the 
3 a 

cr. values compare very well, with an error no greater than 1%. This is most 

likely due to the fact that the model being very "long-coupled" experienced 

negligible pitching; and further, the high cable angle allowed little motion in the 
z direction compared with motion in the x direction. So, the aerodynamic forcing 
effects, which were derived from theory, were small compared with the cable 
tension, which was directly measured. Thus the frequency, a., was largely 
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determined by T q and m, which gave, consequently, less error. Note also 
that for the motion described, aerodynamic damping is largely due to the body* s 
drag: so it would therefore be expected that for a given U, damping would 
increase with a decrease in R . This is, in fact, borne out by the comparison 

3 . 

of the quantitative experimental damping with This type of motion therefore 
has an analogy in airplane long mode "phugoid” oscillations. 

Finally, the theory qualitatively bears out the fact that the system is 
laterally unstable at L = 100 ft. and that it was observed to have a low lateral 
frequency. 

4. 5 Conclusion 

Within the limits of the assumptions listed at the beginnings of Chapters 1 
and 2, the present theory provides a method for predicting the first order motion 
of a large variety of cable-body systems. The key assumptions are that the cable's 
curve must be shallow and that its tension must be essentially constant along its 
length. Otherwise, there are no restrictions on the cable' s motion, i. e. , no 
"instantaneous equilibrium" physical model. Thus the theory may be as readily 
applied to a high frequency system, such as a towed cone in hypersonic flow, as 
to low frequency systems such as towed balloons or the present experimental 
model. 

The essential feature of the theory is that the cable-body system is treated 
as a cable problem, with the body providing end and auxiliary conditions. This 
physical model can lend itself readily to a variety of further applications. For 
instance, the problem of two bodies connected by a cable may be treated by re- 
placing the fixed end condition at s = 0 with a set of end conditions similar to 
those at s = 1, only that these conditions would pertain to the body at that end. 

Another variation could allow for varying cable shape and tension by 
assuming the cable to be composed of finite cable segments each with a given 

A ^ 

T q and a. The equations for each segment are then matched, one to the other, 
through the end conditions of displacement and slope; whereas the end condition 
of the final segment is given by the body equations as before. Similarly, 
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a further application would be to consider a finite body midway along the cable. 

In this case, the end conditions on the two adjacent cable segments are found from 
the equations of motion of the midcable body. 
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APPENDIX I 


The Relation of the Stability Derivatives to the Stability 
Derivatives of Standard Aircraft Convention 


Due to the fact that most information available on stability derivatives 
is based on the x^.y.z^ "wind axis" coordinate system, (Fig. 14) it is 
profitable to give the relations between these and the stability derivatives defined 
in Chapter 2, which are based on the n , n ,n coordinate system. These are 

c x ■ < c x > c V< c z > s \ - t< c x ' + < c z > I S<, „ C V 

u u a a u 

o o o o 

C x . “ (c/bmCjJ c\ + (C Z _) s\ - ((C ) + (C Z ) ]S» o C« o ], 

u u a a u 

o o o o 

c x " < c x > c \ - (C z > s \ + [< c x > - < c z > I s9 o c V 

w a u u a 

o o oo 

C x. - < C/b >« C X.> C \ - < C Z.> ^ e 0 +1<C X.> - < C Z.) ) S<, „ CS o'' 
w “o u o u o “o 

c x -(c/b)[(C z ) C0 O +(C X ) se 0 ], 


C Y = (c ) , 

v p o 


c y. - (C ) , 
V p 0 


C Y =_(C Y ) C0 o +(C Y ) S V 
d d r 


C Y “ -< C Y > C °o - <C Y > S9 o’ 
r r o P o 


c z = < c z » c \ - < c x > s \ - >< c z > - < c x > > S(, o C V 

u u a a u 

o o o o 
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c z . =< c/b >l< c z.> c \-< c x.> s \*« c x.' - (C z.' !S W- 

u u a U a 

o o o o 


c z - <c z > c V (c x > S V [(C Z > + ‘ C X > |S0 O C V 


w a 


u 


u a 

o o 


c .(c/b)i(C ) c 2 « +(C X ) s 2 e o + [(c z _) + (C X .) 1S* 0 C* 0 1. 

* 4 „ 
c, = -<c, ) ce +(C ) so , 

K V ° V ° 

*0 o 


C i. *-<V Ce o*( C n.» S V 
V P P „ 

o o 


C t - < C * > C V< C „ > S \ - » C n » + (C * > ' S8 0 C8 0- 


P P 


P r 

y o O 


C t - < C t » C \ - < C „ > S \ +[<C i ) - V 1S9 ° C V 

r r o P o P o r o 

c =-(c/b)[(c ) ce +(C ) so 1, 

m m o m o 

u u a 

o o 


c = -(c/b) [(c ) ce +<c ) se j, 

m m u o m a 0 

o O 

c = - ( c/b)[(c ) ce + <c ) se j, 

m w m a ° m u ° 

O O 

C m. = Wl( c m .) C0 o+ (C m _) « # |. 

w “o u o 


C m - < o/b > (C m > ' 
q Q o 

C „ - -< C n > C8 o ' <C / > S6 o' 

V P P ^ 

r o o 

C n. * -<V C8 o - <%> S V 

v P P ^ 

r o o 
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c » = <C n > c \ - < c t > s \ + [< c .e > - < c „ 
p p r p 

o o F o 


and 


C n - (C n > c 2 e o +(C , S 2 ^ +[(C ) +<c 

r r o Po P„ 


) JSfl CO , 

r ° ° 

o 

) jse co . 

O 0 

o 
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Fig. 3 Coordinates for the displaced cable. 



Fig, 4 A perturbation of the cable element in the y, £ plane. 



Fig. 5 The body' s coordinate systems. 
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Fig. 8 Coordinate transformation. 



Fig. 9 Sample lateral roots locus plot. 
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Roots Locus Plot for the Lateral 

Case; U=28.0, R =1.09, L = 2.00 
a 



Fig. 18 A lateral roots locus plot. 
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Roots Locus Plot for the Longitudinal 
Case; U = 25. 0, R =1.24, L = 1.76 

cL 



Fig. 19 A longitudinal roots locus plot. 
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tr I o n>|ct 


Tethered Lifting Body Model 



tail airfoil: Clark-Y 

i = 13 1/8 in. f l = 11 in. , i = 8 7/8 in. 

i. Z u 

Fig. 20 Geometry of the test body 



! 



Fig. 21 The test system In the Stanford wind tunnel. 






and 



Q — Rg= 1.24 ft. 
O — R a = 1.09 ft. 
O-R 0 s948 ft. 



Fig. 25 Experimental values of cable critical length, L , vs. wind speed, U. 

cr 


U, (ft /sc 


Fig. 26 Experimental lateral frequency 
wind speed, U. 


□ — R a s 1.24 ft. 
G— R a = 1.09 ft. 
O — R a = 948 ft. 



O, 


o 



BODY ANGLE, a b , (deg) 
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